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I.  STATEMENT  OF  THE  PROBLEM 

"The  general  objectives  of  this  research  program  have  been  the  determina¬ 
tion  and  prediction,  computationally,  of  the  reactive  behavior  of  high-energy 
molecules  and  related  species,  such  as  Intermediates  in  syntheses.  This  has 
Involved  introducing  certain  new  concepts  and  techniques,  as  well  as  extending 
the  range  of  applicability  of  older,  established  procedures.  For  Instance, 
one  of  our  aims  was  to  develop  a  means  whereby  molecular  electrostatic  poten¬ 
tials,  a  very  effective  tool  for  elucidating  electrophilic  reactions,  could 
also  be  applied  to  nrcleophllic  processes. 

He  have  sought  to  achieve  a  better  understanding  of  the  basic  factors 
that  determine  the  reactive  properties  of  certain  important  classes  of  high- 
energy  molecules  (including  nltroaromatlcs  and  nitramines)  aud  how  these  would 
be  altered  by  various  possible  modifications  of  these  molecules.  Our  efforts 
have  also  been  directed  toward  helping  to  design  synthetic  routes  to  other 
potentially-significant  classes  of  energetic  compounds.  Including  the  nitro 
derivatives  of  strained  systems  such  as  cubane. 

-'One  of  our  specific  aims  throughout  this  work  has  been  the  development 
of  a  greater  degree  of  Insight  into  the  factors  governing  the  stabilities  of 
the  compounds  of  interest,  and  how  these  stabilities  are  affected  by,  for 
example,  various  substituents  on  the  molecules.  Our  interests  have  been  in 
both  impact  sensitivity  and  also  in  chemical  stability,  as  toward  moisture. 

'  c4'('cJ?or  U.  ' r  O  L.\ 
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II.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

A.  Development  of  Method  for  Using  Electrostatic  Potentials  to  Predict  and 

Elucidate  Nucleophilic  Reactions 

In  the  past,  the  electrostatic  potential  could  not  be  applied  as  readily 
to  the  analysis  of  nucleophilic  processes  as  to  electrophilic  ones.  This  is 
because  a  positive  potential  does  not  necessarily  indicate  a  corresponding 
tendency  to  react  with  nucleophiles.  The  positive  charges  of  atomic  nuclei, 
being  very  highly  concentrated,  create  strong  positive  potentials  that  may 
outweigh  the  negative  contributions  of  the  electrons,  which  are  dispersed, 
but  nevertheless  may  not  reflect  an  affinity  for  nucleophiles. 
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One  of  our  first  objectives  in  this  project- was  to  find  some  means  for 
overcoming  this  problem,  so  that  nucleophilic  reactions  could  also  be  sub¬ 
jected  to  this  very  useful  mode  of  analysis.  Our  initial  approach  involved 
computing  the  difference  between  the  total  molecular  electrostatic  potential 
and  the  sum  of  the  potentials  of  the  individual  atoms  (taken  to  be  in  their 
free  states)  that  comprise  the  molecule.  It  was  hoped  that  taking  this  dif¬ 
ference  would  eliminate  the  misleading  contributions  to  the  positive  regions. 
This  procedure  was  tested  for  several  molecules  and  found  to  be  partially 
successful  in  predicting  the  sites  for  nucleophilic  attack;  however  it  was  not 
sufficiently  consistent  and  conclusive  to  be  acceptable. 

Accordingly  we  investigated  another  possible  approach,  an  adaptation  of 
a  technique  that  we  Introduced  some  years  ago  for  certain  types  of  electro¬ 
philic  processes  [1-3] .  Rather  than  looking  at  Che  potential  of  the  ground 
state  molecule,  we  allow  it  to  assume  a  configuration  corresponding  to  some 
early  intermediate  stage  in  its  reaction  with  a  nucleophile.  For  an  aromatic 
molecule,  this  means  rotating  the  atom  to  be  replaced  out  of  the  molecular 
plane  and  examining  the  potential  in  the  quasi-tetrahedral  fourth  direction 
that  is  thereby  created.  For  an  aliphatic  system,  this  involves,  for  example, 
computing  the  potential  of  the  structure  in  an  early  phase  of  inversion  of 
configuration.  The  key  to  this  procedure  is  comparing  the  potentials  thus 
obtained  for  the  various  sites  or  systems  of  Interest,  to  determine  which  would 
be  preferred  by  a  nucleophile. 

We  carried  out  several  tests  of  this  idea,  and  found  it  to  work  exactly  as 
desired  [A] .  For  the  fluorobenzene,  j>-nitrofluorobenzene  pair,  for  instance, 
when  each  fluorine  is  rotated  out  of  the  ring  plane,  a  channel  of  positive 
electrostatic  potential  is  created,  a  direct  pathway  whereby  a  nucleophile  can 
approach  the  fluorine  position.  This  pathway  is  considerably  more  positive  in 
the  case  of  p-nitro fluorobenzene,  in  agreement  with  the  experimental  observa¬ 
tion  that  the  presence  of  an  -NO^  group  para  to  a  fluorine  greatly  increases 
the  reactivity  of  the  latter  toward  nucleophilic  substitution.  In  another 
test,  we  applied  this  approach  to  acrolein,  and  were  able  to  correctly  predict 
that  the  carbonyl  carbon  Is  the  most  likely  site  for  nucleophilic  attack. 
Finally,  we  showed  that  methyl  chloride  is  more  susceptible  to  an  S^2  nucleo¬ 
philic  substitution  than  is  methanol,  again  in  full  agreement  with  experi¬ 
mental  results. 
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Thus  Che  proposed  procedure  has  been  successful  In  all  of  Che  cases  inves- 
clgaced.  Ic  of fers  the  possibility,  therefore,  of  opening  a  whole  class  of  | 

important  reactions  to  physically  meaningful  Interpretation  and  prediction. 

Resulting  publication: 

P.  Politzer,  S.  J.  Landry  and  T.  Uamheim,  "Proposed  Procedure  for  Using 
Electrostatic  Potentials  to  Predict  and  Interpret  Nucleophilic  Processes," 

Journal  of  Physical  Chemistry,  86.  4767  (1982). 

B.  Studies  of  the  Electrostatic  Potentials  of  Strained  Systems 

In  view  of  the  Army's  very  considerable  interest  in  nitro  derivatives  of 
strained  molecules  as  potentially  useful  high-energy  systems,  we  have  carried 
out  extensive  and  detailed  studies  of  the  following:  cyclopropane,  triprismane, 
cubane,  homocubane,  j>-blshomocubane,  [l.l.ljpropellane,  and  various  mono-  and 
polysubstituced  -Nf^,  -NH2  end  -CF^  derivatives  of  these  hydrocarbons.  (The 
structures  of  many  of  the  high-energy  molecules  to  be  discussed  are  shown  on 
page  4.) 

A  very  significant  feature  that  we  have  found  to  be  common  to  almost  all 
of  these  strained  hydrocarbons  (one  exception  shall  be  discussed  later)  is  a 
relatively  strong  negative  electrostatic  potential  associated  with  the  C-C 
bonds,  the  most  negative  point  being  near  the  midpoint  of  the  bond.  These  indi¬ 
cate  that  the  bonds  in  these  molecules  can  serve  as  sites  for  electrophilic 
attack,  which  has  indeed  been  observed  in  those  experimental  studies  that  have 
been  made  on  these  systems. 

We  applied  our  technique  of  rotating  hydrogens  away  from  their  equilibrium 
positions  and  comparing  the  strengths  of  the  resulting  electrostatic  potentials 
in  the  open  channels  of  approach  to  the  carbons  that  are  thereby  created. 

(Such  a  carbon  can  be  regarded  as  being  in  an  early  stage  of  a  substitution 
reaction.)  We  Introduced  this  method  for  predicting  the  relative  reactivities 
of  various  sites  about  ten  years  ago  [2,3],  .and  it  has  sines  been  used  success¬ 
fully  by  us  and  by  other  research  groups  [1,4-7].  (In  making  our  predictions,  j 

we  do  also  take  into  account  the  relative  energies  required  for  the  different  1 

hydrogen  rotations.)  In  this  manner,  we  were  able  to  predict  the  most  reactive  ; 

bonds  in  homocubane  [7]  and  in  several  isomeric  blshcmocubanes.  The  former  j 

prediction  is  fully  confirmed  by  laboratory  studies  of  rearrangement  reactions  { 


catalyzed  by  metal  ions.  For  the  bishomocubanes*  relative  bond  reactivities 
have  not  yet  been  determined  experimentally;  however  we  have  presented  our 
findings  and  predictions  at  several  meetings  sponsored  by  Defense  Department 
agencies,  in  order  to  make  them  known  to  researchers  working  with  these  com¬ 
pounds.- 

The  introduction  of  even  one  -N02  group  has  the  consistent  effect  of 
eliminating  the  negative  potentials  associated  with  the  C-C  bonds  in  these 
strained  hydrocarbons;  this  reflects  the  electron-withdrawing  power  of  the 
nitro  group.  With  amine  substituents,  the  situation  is  similar  but  not  as 
extreme;  the  negative  C-C  potentials  are  weakened  but  not  eliminated  18].  The 
configuration  of  the  -NH2  group  in  the  ground  states  of  these  molecules  is 
pyramidal,  'S*N — It  is  interesting  to  note  that  when  this  group  is  forced 
to  be  planar,  —  it  has  the  opposite  effect  upon  the  C-C  potentials;  they 

become  more  negative  than  in  the  parent  hydrocarbon.  This  can  be  interpreted 
a3  involving  a  n-type  charge  donation  by  the  -NH2  to  the  carbon  framework. 

The  preceding  observations  show  the  importance  of  using  correct  geome¬ 
tries,  especially  for  amine  groups.  (Further  examples  will  be  given  later.) 
Accordingly,  we  always  carry  out  an  extensive  computational  optimization  of  the 
geometry  of  each  molecule  that  we  study.  In  a  number  of  cases,  such  as  £-bis- 
horoocubane  and  the  various  aminocubanes ,  there  are,  to  our  knowledge,  no 
experimentally-determined  structures.  Thus  we  provide,  by  our  optimizations, 
the  only  available  structural  data  for  such  systems. 

The  electrostatic  potentials  of  these  substituted  systems  are  also  important 
because  of  what  they  reveal  about  the  substituent  groups  themselves.  We  find 
that  as  the  number  of  -N02  and  -CF3  groups  on  a  given  molecule  Increases,  the 
strong  negative  potential  associated  with  the  oxygens  in  the  former  and  the 
fluorines  in  the  latter  becomes  significantly  weaker.  This  is  important  rela¬ 
tive  to  the  stabilities,  crystal  packing  and  densities  of  the  desired  poly- 
substituted  compounds  in  the  solid  state;  all  of  these  properties  naturally 
depend  upon  repulsion  between  the  polar  substituent  groups.  The  relevance  of 
-CF3  is  that  octaperfluoromethylcubane  has  been  produced  and  studied  in  the 
solid  state.  Thus  our  comparison  of  -N02  and  -CF3  substitution  will  help  in 
anticipating  the  possible  stabilities  and  densities  of  polynitrocubanes,  includ¬ 
ing  octanitrocubane. 


A  strained  molecule  that  has  recently  attracted  a  great  deal  of  interest , 
and  which  shows  some  differences  from  the  others  that  have  been  discussed,  is 
[l.l.ljpropellane.  This  hydrocarbon  is  considered  to  have  a  very  high  degree 
of  strain,  and  it  has  the  particularly  notable  feature  that  all  of  the  bonds  of 
the  bridgehead  carbons  are  in  one  hemisphere;  these  are  the  so-called  “inverted" 
carbons  (see  structure  I) .  Unlike  the  other  strained  molecules  that  we  have 

I 

(inverted  carbons  are 
indicated  by  arrows) 

studied,  the  most  negative  electrostatic  potentials  in  [l.l.ljpropellane  are 
not  associated  with  the  C-C  bonds  but  rather  with  the  bridgehead  carbons.  As 
is  shown  in  Figure  1 ,  there  is  a  relatively  strong  negative  region  to  the  out¬ 
side  of  each  of  these  inverted  carbons;  by  computing  density  difference  plots* 
we  confirmed  that  this  is  due  to  a  buildup  of  electronic  charge  in  each  of  these 
regions,  which  is  accompanied  by  a  decrease  (relative  to  the  superposed  free 
atoms)  between  the  bridgehead  carbons.  Thus  we  predict  that  the  most  reactive 
sites  toward  electrophiles  should  be  to  the  outside  of  each  bridgehead  carbon. 

The  available  experimental  data  concerning  the  reactive  behavior  of  [1.1.1]- 
propellane  are  in  accord  with  this  prediction  (9).  We  also  anticipate  some 
degree  of  diradical  character,  again  focusing  upon  these  regions  by  the  inverted 
carbons. 

Resulting  publications: 

(1)  P.  Politzer,  L.  N.  Doraelsmith,  P.  Sjoberg  and  J.  Alster,  "Electrostatic 
Potentials  of  Strained  Systems:  Nitrocyclopropane",  Chemical  Physics  Letters, 

92,  366  (1982). 

(2)  P.  Politzer,  L.  N.  Domelsmith  and  L.  Abrahmsen,  "Electrostatic  Poten¬ 
tials  of  Strained  Systems:  Cubane,  Homocubane  and  Bishomocubane",  Journal  of 
Physical  Chemistry,  88,  1752  (1984). 

(3)  P.  Politzer,  K.  Jayasuriya  and  B.  A.  Zilles,  "Some  Effectsof  Amine 
Substituents  in  Strained  Hydrocarbons",  Journal  of  the  American  Chemical  Society, 
in  press  (1985) . 


Figure  X,  Calculated  electrostatic  potential  of  [l.l.ljpropellane, 

as  computed  from  an  ab  initio  SCF  ST0-5G  wave  function 

* 

using  a  geometry  that  had  been  optimized  at  the  6-31G 
level  with  the  GAUSSIAN  80  program*  The  potential  is  shbwn 
in  the  plane  passing  through  the  two  bridgehead  (Inverted) 
carbons  and  the  methylene  carbon.  Contour  values  are  in 
kcal/mole.  The  positions  of  the  most  negative  values  of 
the  electrostatic  potential  (potential  minima)  are  indicated 
by  the  solid  circles,  and  correspond  to  -30.0  kcal/mole. 


Introduction  and  Applications  of  the  Bond  Deviation  Index 


Early  in  this  project,  we  began  to  use  the  concept  of  "bond  path",  pioneered 
by  Bader  and  his  colleagues  [10],  as  a  means  of  studying  chemical  bonds.  The 
bond  path  is  the  ridge  of  maximum  electronic  density  that  links  two  chemically-* 
bonded  nuclei.  (In  cyclopropane,  for  instance,  the  C-C  bond  paths  curve  out¬ 
ward,  consistent  with  the  idea  of  bent  bonds;  in  propane,  on  the  other  hand, 
they  are  straight  [10,11].)  In  order  to  put  this  on  a  more  quantitative  basis, 
we  have  introduced  a  new  property,  the  "bond  deviation  index"  [12].  This  is  a 
numerical  measure  of  the  extent  to  which  an  actual  bond  path  deviates  from  a 
reference  path  defined  in  terms  of  the  superposed  electronic  densities  of  the 
undistorted  free  atoms  placed  at  the  same  positions  as  they  occupy  in  the  mole¬ 
cule.  As  examples,  the  bond  deviation  indices  for  the  C-C  bonds  in  propane  and 
cyclohexane  are  0.000;  in  cyclopropane,  they  are  0.080,  and  in  tetrahedrane, 
0.120  [12]. 

While  many  of  our  applications  of  this  type  of  analysis  have  been  to 
strained  molecules,  it  is  by  no  means  limited  to  such  systems,  and  we  have  also 
used  it,  for  instance,  to  study  the  bonding  in  certain  organophosphorus  com¬ 
pounds.  The  bond  deviation  index  represents  a  quantitative  means  of  character¬ 
izing  chemical  bonds  and  allowing  meaningful  comparisons  to  be  made.  For 
example,  our  calculated  bond  deviation  indices  for  triprismane  show  that  its 
C-C  bonds  can  be  viewed  as  being  of  two  types:  cyclopropane-like  (the  bonds 
in  the  three-sided  faces),  and  cubane-like  (the  bonds  connecting  the  three- 
sided  faces) .  The  calculated  electrostatic  potentials  of  the  triprismane  C-C 
bonds  confirmed  the  validity  of  their  being  classified  as  belonging  to  these 
two  categories  [8]. 

A  key  advantageous  feature  of  the  bond  deviation  index  is  that  it  permits 
one  to  focus  specifically  upon  any  individual  bond  in  a  molecule,  rather  than 
being  obliged  to  deal  with  quantities  that  reflect  the  molecule  as  a  whole.  By 
means  of  the  bond  deviation  index,  one  can  observe  and  compare  the  effects  of 
chemical  substitution  and  structural  distortion  upon  each  individual  bond  in  a 
molecule.  We  have  now  computed  these  indices  for  the  bonds  in  about  30  mole¬ 
cules,  and  have  used  them  in  analyzing  the  reactive  properties  of  these  systems. 
For  example,  they  were  an  integral  part  of  the  studies  (mentioned  above)  of  the 
reactive  behavior  of  cubane,  honocubane  and  £-bishomocubane  [7]. 


There  are  also  Indications  that  the  bond  deviation  index  ma.,  provide  a 
means  for  predicting  when  the  total  strain  energy  of  a  molecule  is  the  sum  of 
the  strain  energies  of  its  component  parts.  For  instance,  as  mentioned  above, 
the  Indices  for  trlprismane  show  that  it  has  two  cyclopropane-like  and  three 
cubane-llke  faces.  And  Indeed,  taking  twice  the  strain  energy  of  cyclopropane 
plus  3/6  times  that  of  cubane  (  which  has  six  faces)  actually  gives  a  result 
that  is  close  to  the  strain  energy  of  trfprismane  [13]. 

Resulting  publication: 

F.  Politzer,  L.  Abrahmsen,  F.  Sjoberg  and  F.  R.  Laurence,  "The  Use  of 
Bond  Paths  for  Quantitatively  Characterizing  Bond  Strain",  Chemical  Physics 
Letters,  102,  74  (1983). 

Results  and  discussions  of  bond  deviation  index  calculations  have  also  been 
Included  in  several  others  of  our  publications,  listed  elsewhere  in  this  report. 

D.  Studies  of  Mltroaromatlc  Systems 

We  have  computed  wave  functions,  structures  and  various  properties  (electro¬ 
static  potentials,  electronic  densities,  force  constants,  energies,  atomic 
charges,  etc.)  for  about  15  mono-  and  disubstltuted  benzene  derivatives.  Most 
of  these  included  at  least  one  -N02  substituent,  with  -NH2,  -CHj  and  -NO  as  other 
possibilities.  We  have  used  these  results  to  study  the  effects  that  these  sub¬ 
stituents  have  upon  each  other  and  upon  the  overall  reactivities  and  stabilities 
of  the  molecules. 

For  Instance,  by  comparing  our  computed  C-N02  stretching  force  constants 
for  nitrobenzene,  m-dlnltrobenzene,  and  o-  and  p-nltroanlline,  we  have  ascer¬ 
tained  that  the  C-NO^  bond  is  strengthened  by  the  presence  of  an  -NH2,  parti¬ 
cularly  when  the  latter  is  in  the  ortho  position,  but  is  relatively  unaffected 
by  a  meta-NOy.  We  have  also  investigated  how  rotation  of  the  -N02  in  nitro¬ 
benzene  around  the  C-N  bond  affects  the  reactivities  of  various  sites  on  the 
ring.  While  the  most  stable  orientation  for  the  nitro  group  in  this  molecule 
is  in  the  plane  of  the  ring,  steric  factors  may  force  it  to  rotate  out  of  this 
plane  in  some  polysubstituted  systems.  This  affects  its  conjugation  with  the 
ring,  and  may  change  the  reactivities  of  various  sites.  Out  preliminary  conclu¬ 
sion  is  that  this  does  occur,  but  that  the  effect  is  a  relatively  small  one. 

The  influence  of  the  nitro  group  is  exerted  primarily  inductively,  and  this  does 
not  involve  conjugation  with  the  ring. 
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The  strong  electron-withdrawing  power  of  the  -NO^  group,  which  has  manifest¬ 
ed  Itself  by  eliminating  the  negative  electrostatic  potentials  associated  with 
the  C-C  bonds  in  strained  hydrocarbons  (see  above),  is  now  found  to  eliminate 
the  negative  regions  above  and  below  the  aromatic  ring  that  are  present  in  ben¬ 
zene  and  especially  in  aniline,  reflecting  their  it  electrons  [14].  In  the  case 
of  the  nltroanillnes,  which  we  have  examined  in  particular  detail,  we  have  found 
that  the  power  of  the  -NO2  group  to  deactivate  the  aromatic  ring  toward  electro¬ 
attack  overwhelms  the  opposing,  activating  tendency  of  the  -NH2  [14].  To  the 
extent  that  the  nltroanillnes  do  undergo  electrophilic  substitution,  however, 
our  potential  plots  Indicate  that  it  should  be  governed  by  the  directing  proper¬ 
ties  of  the  amine  groups. 

These  studies  again  underscore  the  Importance  of  geometry  optimization.  In 
aniline,  we  find  the  -NH^  group  to  be  pyramidal,  in  agreement  with  experimental 
determinations;  in  the  nltroanillnes,  however,  our  optimizations  show  it  to  be 
planar.  The  strong  electron-withdrawing  effect  of  -N02  apparently  causes  the 
-NH2  to  adopt  the  planar  configuration,  in  which  it  is  better  able  to  act  as  a 

n-electron  donor,  via  Increased  contributions  from  structures  such  as  II, 

2 

in  which  the  amine  nitrogen  is  sp  -hybridized,  with  its  hydrogens  in  the  plane 
of  the  ring. 

H  <+)  /=\  (+)  0,"> 

11  h'N 


Resulting  publication: 

P.  Politzer,  L.  Abrahmsen  and  P.  Sjoberg,  "Effects  of  Amino  and  Nitro 
Substituents  upon  the  Electrostatic  Potential  of  an  Aromatic  Ring",  Journal  of 
the  American  Chemical  Society,  106,  855  (1984). 

E.  Discovery  of  Positive  Electrostatic  Potential  Buildups  Near  Aromatic 


C-NO^  Bonds 

In  the  course  of  the  studies  of  aromatic  systems  described  above,  we  have 
discovered  that  the  electrostatic  potentials  of  nitroaromatic  molecules  show 
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unexpected  buildups  of  positive  potential  above  and  below  the  C-N02  bond  regions. 
Their  magnitudes  depend  upon  the  natures  and  locations  of  other  substituents  on 
Lhe  aromatic  ring.  Such  effects,  when  not  at  nuclear  positions,  are  certainly 
uncommon  and  we  have  accordingly  investigated  this  in  some  detail  in  the  case 
of  nitrobenzene.  On  the  basis  of  its  computed  electric  field  and  also  calcula¬ 
tions  of  its  interaction  with  OH  ,  we  have  reached  the  very  intriguing  conclu¬ 
sion  that  the  C-N02  bond  in  nitrobenzene  (and  presumably  in  the  other  nitro- 
aromatics  as  well)  can  serve  as  an  initial  site  for  nucleophilic  attack. 

Showing  the  key  role  played  by  electrostatic  effects  in  the  OH”  interaction  is 
the  fact  that  we  find  a  direct  relationship  between  the  potential  of  nitroben¬ 
zene  at  various  points  in  the  C-NO2  region  and  its  Interaction  energy  with  OH** 
at  those  points. 

Some  immediate  mechanistic  applications  of  these  results  are  in  the  reac¬ 
tions  of  nltroaromatic  molecules  with  hydroxide  and  alkoxlde  ions  to  form 
Meisenheimer  complexes  [15],  and  in  their  reactions  with  hydrogen  halides  that 
are  being  studied  by  Nielsen  and  co-workers  [16] .  The  ability  of  C-N02  bond 
regions  to  serve  as  initial  sites  for  nucleophilic  attack  may  well  be  a  very 
significant  contributing  factor  to  whatever  degrees  of  chemical  instability 
are  shown  by  various  nltroaromatic  molecules,  as  toward  hydrolysis  for  instance. 

Resulting  publication: 

P.  Politzer,  P.  R.  Laurence,  L.  Abrahmsen,  B.  A.  Zilles  and  P.  Sjoberg, 

"The  Aromatic  C-NC^  Bond  as  a  Site  for  Nucleophilic  Attack",  Chemical  Physics 
Letters,  111,  75  (1984). 

F.  Mechanistic  Studies  of  Aromatic  Nitration 


In  order  to  help  clarify  some  aspects  of  the  mechanisms  of  aromatic  nitra¬ 
tions,  we  have  computed  the  structures,  stabilization  energies,  and  atonic 
charge  redistributions  associated  with  the  formation  of  both  weakly-  and 
strongly-bound  benzene  «N0^  and  toluene* NO ^  complexes  that  are  possiole  Inter¬ 
mediate  stages  in  the  nitration  of  these  aromatic  molecules.  The  Initial  weakly 
bound  complexes  are  relatively  delocalized,  with  the  NO^  approximately  parallel 
to  the  aromatic  ring  and  roughly  2.8  A  above  it,  in  a  region  of  negative  electro 
static  potential.  The  strongly-bound  complexes  are  localized  at  specific  ring 
carbons  (three  Isomers  in  the  case  of  toluene),  and  can  be  regarded  as  sigma 
complexes . 


Certain  Interesting  and  consistent  patterns  can  be  seen  in  our  results. 

One  of  these  is  the  fact  that  the  same  type  of  resonance  structure  appears  to  be 
the  dominant  contributor  to  each  of  the  four  sigma  complexes.  Another  is  the 
picture  that  emerges  of  the  two  functions  performed  by  the  methyl  group  in 
toluene  and  its  complexes:  First,  it  induces  a  rearrangement  of  charge  within 
the  aromatic  ring  whereby  the  carbon  to  which  it  is  attached  becomes  more  posi¬ 
tive  (in  agreement  with  the  NMR  spectra  of  benzene  and  toluene),  by  essen¬ 
tially  the  same  increment  in  each  case.  The  resulting  electrostatic  interaction 
C+,^  may  be  responsible  for  the  various  stabilization  effects  in  the 
toluene -NO^  sigma  complexes  that  result  in  the  well-known  ortho,  para  directing 
properties  of  -CHj.  A  second  function  of  the  methyl  group,  under  the  stimulus 
of  the  strong  electron-attracting  power  of  the  -NC^,  is  as  a  source  of  elec¬ 
tronic  charge,  which  comes  from  its  hydrogens.  This  is  observed  in  the  weakly- 
bound  complex  as  well  as  the  strongly-bound  ones. 

Resulting  publication: 

P.  Politzer,  K.  Jayasuriya,  P.  Sjoberg  and  P.  R.  Laurence,  "Properties  of 
Some  Possible  Intermediate  Stages  in  the  Nitration  of  Benzene  and  Toluene", 
Journal  of  the  American  Chemical  Society,  in  press  (1985). 


G.  Relationships  for  Correlating  and  Predicting  Impact  Sensitivities  of 
Nltroaromatic  Compounds 

He  have  discovered  two  good  correlations  between  the  experimentally-measured 
impact  sensitivities  of  a  group  of  nltroaromatic  compounds  and  certain  electronic 
properties.  The  compounds  are:  2,3,4,6-tetranitroaniline,  2,4,6-trinitrophenol, 

1.3.5- trinitrobenzene,  2,4 ,6-trinitroaniline,  2,4,6-trinitrotoluene,  2,4-diamino- 

1.3.5- trinitrobeozene  and  2,4, 6- triamino-1, 3 ,5- trinitrobengene. 

First,  there  is  a  direct  relationship  between  the  impact  sensitivity  and 
the  positive  potential  at  the  C-NC>2  bond  midpoint,  Vn^,  that  can  be  computed 
from  the  calculated  atomic  charges  of  the  carbon  and  nitrogen  atoms  [17] .  (Both 
of  these  charges  are  generally  positive.)  The  relationship  is  shown  in  Figure  2. 

Second,  there  is  a  good  correlation  between  impact  sensitivities  and  the 
unbalanced  electrostatic  forces  within  the  -NO2  groups  in  the  crystalline  forms 
of  the  compounds  (Figure  3) .  The  unbalanced  forces  come  about  because  the 
crystalline  environments  distort  the  molecules  from  their  Isolated  gas-phase 
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equilibrium  geometries.  As  a  result,  the  net  intramolecular  electrostatic  force 
upon  each  nucleus  (due  to  the  electrons  and  other  nuclei  in  the  molecule)  is  not 
zero,  as  it  should  be  at  equilibrium,  but  rather  has  some  finite  value.  We 
calculate  these  forces  from  the  molecular  wave  function  that  is  computed  using 
the  molecule's  geometry  in  the  crystal.  By  summing  the  net  forces  that  are 
exerted  upon  all  of  the  -NO2  oxygens  in  a  molecule,  we  obtain  a  quantity  that 
correlates  well  with  the  compound's  measured  Impact  sensitivity.  Equally  good 
results  are  produced  by  including  the  net  forces  upon  the  -NO^  nitrogens  in  the 
summations  (Figure  3) . 

These  relationships  provide  insight  into  some  of  the  factors  that  determine 
impact  sensitivities  and  what  consequences  may  be  anticipated  from  modifications 
such  as  different  functional  groups.  They  will  also  allow  reliable  predictions 
to  be  made  for  compounds  that  have  not  been  tested;  indeed,  the  correlation 
involving  Vm^n  can  be  used  in  evaluating  potentially  important  compounds  that 
have  not  yet  even  been  synthesized,  since  it  does  not  require  the  crystalline 
structure. 

H.  Studies  of  Hltramlnes 

In  the  course  of  a  detailed  study  of  nitramines,  we  have  computed  optimized 
structures,  energies,  electrostatic  potentials,  bond  deviation  indices,  bond 
force  constants,  atomic  charges  and  overlap  populations  for  more  than  a  dozen 
saturated  six-membered-ring  systems,  containing  three  or  fewer  ring  nitrogens 
and  substituted  with  zero,  one,  or  tvo  nltro  groups.  Many  of  these  structures 
were  not  previously  known.  The  interest  in  these  molecules  is  of  course  due  to 
their  relationship  to  RDX  and  HMX,  currently  two  of  our  most  effective  explo¬ 
sives.  For  somewhat  similar  reasons,  we  have  also  computed  the  properties  of 
hexamethylenetetramine  (hexamine) ,  which  is  a  starting  material  in  the  manufac¬ 
ture  of  both  RDX  and  HMX. 

Using  these  results,  we  are  analyzing  the  reactivities  and  stabilities  of 
these  systems,  and  how  they  are  affeated  by  structural  modifications,  such  as 
varying  the  number  of  ring  nitrogens  and  the  number  and  locations  of  -HOj  groups. 
We  have  been  able  to  show,  for  example,  that  the  N-N  bond  is  stronger  in  dimethyl 
nltramlne  than  in  cyclic  nitramines,  which  is  in  agreement  with  recent  findings 
in  laser- induced  decomposition  studies  [18].  One  of  our  alms  in  this  work  is 
to  elucidate  the  decomposition  pathways  of  nitramines,  such  as  RDX  and  HMX. 
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I.  Development  of  Method  for  Calculating  Molecular  Energy  Differences  froa 
Electrostatic  Potentials  at  Nuclei 

In  the  course  of  our  work  with  electrostatic  potentials,  we  developed  a 
formula  that  gives  the  difference  In  energy  of  tvo  isoelectronlc  systems  In 
terms  of  the  electrostatic  potentials  at  their  nuclei.  When  tested  for  nine 
pairs  of  diatomic  molecules  and  molecular  ions,  the  average  error  relative  to 
the  experimental  value  was  only  0.281.  It  is  particularly  noteworthy  that  in 
every  Instance,  our  energy  difference  is  more  accurate  than  the  Hartree-Fock 
prediction,  even  though  our  potentials  were  obtained  from  the  same  Hartree-Fock 
wave  functions.  An  immediate  application  of  our  formula  is  in  the  calculation 
of  molecular  electron  affinities,  which  are  often  quite  difficult  to  determine 
experimentally.  Overall  good  results  were  obtained  in  computing  the  electron 
affinities  of  Cj,  OF  and  CN. 

Resulting  publication: 

P.  Polltzer  and  P.  Sjoberg,  "A  Formula  for  Calculating  Molecular  Energy 
Differences  from  Electrostatic  Potentials  at  Nuclei",  Journal  of  Chemical 
Physics,  78,  7008  (1983). 
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